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Abstract 



Simple bosonic path integral representation for path ordered ex- 
ponent is derived. This representation is used, at first, to obtain new 
. variant of non-Abelian Stokes theorem. Then new pure bosonic world- 

line path integral representations for fermionic determinant and Green 
functions are presented. Finally, applying stationary phase method, 
we get quasiclassical equations of motion in QCD. 



1 Introduction 

Elimination of fermionic degrees of freedom is very desirable in many prob- 
lems of quantum field theory and elementary particle physics. First, in lattice 
gauge theories integration with respect to Grassmannian variables leads to 
serious complications of numerical simulations. The second, appearance of 
fermionic variables in functional integrals hampers the application of station- 
ary phase method. As a consequence, one cannot also apply quasiclassical 
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expansions to evaluation of functional integrals in theories with fermions ex- 
cept some especially simple cases. 

Indeed, typical Green function can be written as 



G{x x ,...,x n ) = J DAD^D^e^ SYM W +s ^*>*' A »0 1 (x 1 )...O n (x n ) (1) 

where Oi are some operators whereas Sym and Sf erm are Yang-Mills and 
fermionic actions respectively Quasiclassical approximation is defined, up to 
some subtleties, by stationary point equations 

-j^- = [the source of YM field] (2) 
But what must be written in R.H.S. of eq. (0) ? Of course, one cannot 

put 

[the source of YM field] = 5S ^ m (3) 

because Sf erm depends on Grassmanian variables. Moreover, without pre- 
liminary exception of fermionic variables one cannot write in R.H.S. of (0) 
nothing except zero. But this means that in zero approximation YM field can 
be considered as free. It seems inappropriate in all cases in which interaction 
is strong. 

So for application of quasiclassical methods as well as for facilitation of 
numerical simulations on the lattice fermionic variables in functional integrals 
of the type ([!]) must be integrated out and result must be represented as 
functional integral with respect to only bosonic variables. In other words, 
the theory must be bosonized. 

The problem of bosonization of fermionic theories has a long history. 
Most likely, the first example of bosonisation of fermionic theory was given 
by Schwinger in his famous paper |T[ concerning full solution of massless 
QED 2 . Then the problem of bosonisation was investigated by many authors, 
but more or less exhaustive solution was obtained only in two dimensional 
case and in some three dimensional models (see, for instance, papers and 
references therein). In realistic four dimensional case only partial success 
was achieved (see, for instance, 0). In fact in all proposed bosonization 
schemes in four dimensions it is necessary to evaluate (exactly or in some 
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approximation) fermionic determinant - but it is just the main problem that 
must be solved by means of bosonization. 

To author's knowledge, the only exceptions are recent papers by Lusher 
[|] and Slavnov || (see also ||). In Ref. [|J fermionic determinant on the 
finite lattice is represented as infinite some of bosonic determinants. In Ref. 
U fermionic determinant in D dimensions is expressed via bosonic one in 
D + 1 dimensions. These approaches seem useful in lattice theories but they 
cannot be applied to investigation of quasiclassical approximation. 

So hitherto no quite satisfactory representation for fermionic determinant 
in terms of bosonic fields is known, and at present paper we will develop an- 
other approach to bosonization. Namely, we will derive pure bosonic world- 
line path integral representation for fermionic determinants, Green functions 
and Wilson loops. 

Worldline approach to quantum field theory also has very long history. 
It was originated many years ago in classical works by Feynman and 
Schwinger || . The main idea of this approach is to represent fermionic 
determinants and fermionic Green functions as functional integral over tra- 
jectory of a single relativistic particle. 

Let us consider, for instance, fermionic determinant for SU(N) Yang- 
Mills theory in Euclidean space: 



where V = 7^V M = 7^ - iAJ, (7^ = 7^, {7^, 7"} = ^\ iA^x) G 
su(N). 




(4) 



One can write: 




-lndet(-V At V M + a^F^ + m 2 ) (5) 



where 




(6) 




(7) 



Further, eq. (|5p can be written, up to inessential constant, as 
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InD = -tr ln(- V^ t + a^F^ + m 2 ) 

= _l r^ e -™\ e -n-v,v^) (8) 

The integral (§) is divergent at T = and must be regularized. One can 
use, for instance, ^-function regularization: 



2 ds T(s) Jo s=o V ' 

However, in what follows we shall write, for short, formal expression (H) for 
In D bearing in mind any suitable regularization. 

Further, trace in eq.([|) can be represented as functional integral: 



= J pBc DqtTPe^S.- £ dt^-i^A,(q)+Ta^F, u (q))\ (10) 

Here PBC means 'periodic boundary conditions'. Path ordering in (|T0|) cor- 
responds to both spin and colour matrix structures. 

Let us suppose for a moment that spin and colour degrees of freedom 
are absent in fllPf ). (This is just the case of scalar electrodynamics). Then, 
integrating in (HD over fermionic fields and using for arising fermionic deter- 
minants and fermionic Green functions formulae of the type ([|), (|10D , one 
obtains formulation of quantum field theory in terms of particles interacting 
with gauge field A. This is just result of classical work bu Feynman 0. 

In realistic models, however, one must take into account spin and colour 
degrees of freedom. So, to develop worldline formulation of quantum field 



theory, it is necessary to represent the path ordered exponent in fliTf ) as func- 
tional integral. In the case of QED this was done by Fradkin in terms 
of fermionic path integral. This representation and their modifications were 
successfully used, in particular, for construction of derivative expansion in 
QED and for investigation of complicated Feynman diagrams [0. It is 



also intensively used for investigation of hidden supersymmetry in fermionic 



4 



theories (see, for example, [|T2[ and references therein). Recently, D'Hoker 
and Gagne derived fermionic path integral representation for fermionic de- 



terminants for particles coupled with arbitrary tensor field [13 . 

In papers fl8f-||13| fermionic path integral representations were derived 
and used only for spin degrees of freedom. For colour P-exponent in (|10|) 
fermionic path integral representation also can be derived (see, for example, 
though it seems not so elegant as one for spin P-exponent. Apropos, 



works |L4| were, to author's knowledge, the first attempts to obtain non- 



perturbative information about QCD in worldline formalism. 

But, as we already stated at the beginning of the present paper, there 
exist important problems for which pure bosonic worldline path integral rep- 
resentation for fermionic determinants and Green functions is very desirable. 
We see that for solution of the latter problem it is sufficient to derive bosonic 
worldline path integral representation for the trace of path ordered exponent 

Z = tiPe i ^ dtB ^ (11) 

with matrix B(t) G GL(N), B(0) = B(l). Indeed, substituting such repre- 
sentation with 

B(t)=q»A,{q(t))-o^F, v {q(t)) (12) 

in (|lOl), one obtains desired representation for fermionic determinant. As 
we will show later, bosonic path integral representation for Z allows also 
to obtain bosonic worldline path integral representation for fermionic Green 
functions. 

Different bosonic path integral representations for Z in the case B(t) G 
SU{2) were proposed by several authors (see P~5|-[]TS||) . For more general 



case B(t) G su(N) the an analogous representation was pointed out in [[19 
The typical result for B G SU (2) is 



DS(t) exp\- / dt 
2 jo 



tr<r 3 (S(t)B(t)S^(t) + iS(t)&(t) 



(13) 



The integration in ([0|) is carried out over all trajectories in the group 
SU(2), S G SU(2), a 3 is a Pauli matrix. 

The last result has two disadvantages. First, in any parametrization of 
SU(2) the "action" in (p~3|) is rather complicated non-polinomial function. 
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This hampers the usage of (0) in practice. The second, it appears that inte- 



gral (13) is ill-defined and needs insertion of some regulators in the "action" 



(see the discussion of this point in |T6|, [0| and in more recent paper [FZD] ) . 

Recently Dyakonov and Petrov found much more elegant formula for path 
ordered exponent. Namely, for the case when Z is Wilson loop and gauge 
group is SU(2), they derived from (|i~3"|) the following expression for Z: 



j-. i £ dx^ Afj,(x) 

tr F e j t 

J Dn{x) [] (n a {x)n a {x) - 1) 



exp ji jf dx» A dx u [-F; v n a + e ahc n a D^n b D u n c ] j (14) 

where £ is two dimensional surface spanned on contour 7, a = 1, 2, 3, - 
covariant derivative. This formula can be considered as non-Abelian variant 
of Stokes theorem. We shall continue the discussion of this result in section 4. 

In the present paper we will derive alternative bosonic path integral rep- 
resentation for Z in the general case B G su(N). The "action" in this 
representation is quadratic and so it is much more simpler then one in (O) 



and fllil) . This derivation is presented in the section 2. In the section 3 
we check the representation obtained by direct evaluation of the functional 
integral that defines Z. In fact, we give alternative proof of results obtained 
in the section 2. In the section 4 we derive an non-Abelian analog of Stokes 
theorem and compare our results with those due to Dyakonov and Petrov. 
In the section 5 we derive bosonic worldline path integral representation for 
fermionic determinant and Green functions in Euclidean QCD. In the section 
6 we get analogous results for Minkowski space and then, applying stationary 
phase method, derive quasiclassical equations of motion in QCD. In the last 
section we summarize our results and discuss perspectives of future investi- 
gations. In two appendixes we derive some auxiliary formulae that are used 
in the main text of the paper. 



2 Bosonic path integral representation for the 
trace of path ordered exponent 

Let N x N matrix U(t), < t < 1, be defined by equations 
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U(0) = I N (15) 
where In is N x N unit matrix. Then, obviously, 

Z = trU(l) (16) 
First, we consider the case B 6 su(N), that is 

B^ = B, trfi = (17) 
Let us consider an operator 

B = a\.B T s a s (18) 

that acts in some Fock space T. Operators a\ and a s in (|T8| ) are usual bosonic 
creation and annihilation ones, [a r ,aj] = 5 r s . 

Let Ti n be n-particle subspace of J 7 , U n is orthogonal projector on Ti n . 
Then 

Z = trn 1 Pe i /o d '^W (19) 
Indeed, if N — aja r , then 



B,N 



(20) 



So -B7i n C 7i n and in one dimensional subspace TC\ the operator B can 
be identified with the matrix B via relation: 

Bat<f> r \0>= al(B r s <f> s )\0 > (21) 
Projector IT can be represented as 



IT = / £ d\5(N~\) 

Jl-e 

rfA [ ^L e -iM+mv (22) 



l-e 



2tt 



-<x> 
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where < e < 1. 
So 



Z = lim trn ie -^Pe^o d *^) 

5^+0 

oo 

= lim [ dX [ ^tre»»«*- s >*e i fi* A ® 

5^+0 Jl-e J 2ir 



-e 

— oo 
oo 



= lim / " dX [ ^tre-^trPe'/o^+^W (23) 

— oo 

The latter equality is valid due to (|20|). 



The trace of ordered exponent in (23) can be represented as functional 



integral. In our case its explicit form essentially depends on some subtleties 
in its definition. So let me remind in a few words the general construction of 



functional integral. For more detailed discussion see, for instance, p2 
Let 

Z = tiVe^ dt& (24) 

where H = H(a' , a; t) is an operator in Fock space T. 
In our case 

H = (irj — 8)N + %B (25) 
Equivalently, we can write 

H = H(p,q;t) (26) 

a + a — , 

" = ^T' * = -S?T (27) 

Let H = H(p, q; t) = H(z, z; t) be Weyl, or normal, or any other symbol 
of the operator H, 

P-iq - P + iq / 0fi N 

2 = ^T' z = ~n (28) 
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Further, let * be the operation that represent the multiplication of op- 
erators in the language of symbols. This means that if K = K(p, q) and 
L = L(p, q) are some operators with symbols K = K(p, q) and L = L(p, q) 
and M = KL then 

M(p,q) = (K*L)(p,q) (29) 

In this terms 

Z = lim / dpdq ( e « H(v;0) * e » H(,,,5 » ) * ... * e^''" ,s ^A (p, q) (30) 

n — *oo J V / 

The last formula can be rewritten as functional integral. Its concrete form 
depends on the choice of a kind of symbols used in (|3~Cf) . In particular, for 
Weyl symbols 

Z = Af[ DpDqefo dt ( ip MiW +Hw< *(t),<im)) ( 31 ) 

JPBC 

whereas for normal symbols 

Z = lim Af' [ Dzm J' ) dt{z(t)z(t)+H n or m (mAt+^)) (32) 

JPBC 

Here PBC means 'periodic boundary conditions', z and z are independent 
complex variables, and A/", Af' are normalization constants (that we will usu- 
ally omit in what follows). 
In our case 

H w (p,q;t) = H w (z\z;t) 

N 

= (irj — 5)z' z + iz< Bz — —(irj — 5) (33) 
where p, q are connected with z by formulae (|28|) , and 

H norm (z, z) = (ir] - 5)z ] z + iz ] Bz (34) 

Formulae (|3l"l) , ( |32"|) correspond to standard sign conventions. However, 
for our purposes it is convenient to change variables 

z{t) -> z(l - t) 
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Then formulae ([H]), (0) with symbols ([33|), (0) can be rewritten as 



D z exp < i / rft 

'PBC I JO 

+ AT 



zz f (t)i(t) + z\t)B(t)z(t) 



where -D 2 2: = .D(Rez).D(rmz) and 



(35) 



lim 



£>z.Dz exp ^ / dt [iz(t)z(t) + z(t)B(t)z(t - e) 



+{r] + iS)z{t)z{t - e)}} (36) 



respectively. 

There are several essential differences between formulae ([35|) and 
First, in (35) z and z' are complex conjugated variables whereas in 



1 
161) 



2; 

and z are independent. The second, the last term in the "action" in (^) is 
absent in (0). The third, there is the shift of "time" variable in the last two 
terms in the "action" in (^) that is absent in (|3~o]). In the next section we 
will show by explicit calculations that this shift just compensates the absence 
of the term 

in (^). Finally, it is worth to note that the limit in (|36|) must be evalu- 
ated after functional integration in ( j36|) because this two operations don't 
commute. It will be confirmed by explicit calculations in the section 3. 

Substituting (|3~5|) and (|36]) in ( p3|) and evaluating limits e — > 0, 5^0 
(bat not the limit e — > +0!), ones obtains, respectively, 



A 

D 2 z J ^l^fodtiizU+z^Bz+nizU-l-^) ^ 



and 



10 



°° d 1 
Z = lim / DzDz / — exp \ i / dtizz + zBe~ e ^ z 

e^+OjpBC J 2n [Jo 

— oo 

+ r/(ze-^z-l)} (38) 

We won't try to justify here the validity of limiting procedure e — > 0, 
5^0 because in the next section we will check formulae (]3~7D , fl3~8f ) by direct 
calculation. 

In ( [37D one can integrate with respect to rj: 

Z = J pBc D 2 z5(£ dtzh - 1 - ^ylo^^Bz) (3g) 
But we cannot obtain 5-function by integration with respect to 7] in (|38|) 



because z and z in (|38|) are independent complex variables. 



One can also get another useful form of the representation ( ^9j), namely 
Z = [ D 2 z l[6(z\t)z(t) - 1 - EySl^z+z^Bz) ( 4Q ) 

J PBC -j. 2 

In what follows, we will sometimes omit symbol T7 m formulae of the type 
©• 

To get formula (pEOP, it is sufficient to insert projectors 111 represented in 



the form (|22| ) between each pair of adjacent factors in (^) and to repeat all 
calculations that have led us to representation (|3"§D. 

The same arguments allow also to obtain an analogous variant of ([38]): 

Z = lim / Dz(t)Dz(t)Dr](t)exp{i f 1 dtiz(t)z(t) + z(t)B(t)e' e ^z(t) 

c^+OJpBC JO 

+ r,{t){z(t)e**z{t)-l)} (41) 

Thus we have derived four different representations (|3"8"|)-(f4"l~D for one 
quantity Z. They all appear to be useful in quantum field theory. 

We have got (^1)-(41) assuming that B 6 su(N). But these representa- 



tion remains valid for every trace free matrix by virtue of analytical contin- 
uation. Representation for B e GL(N) can be obtained by extracting of the 
trace part of the matrix B at the beginning of calculations. Indeed, if 
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B = B' + -I N tvB 



then tr B' = and 



tr Pe l fo dtB = e ^ dt tr B tr P e l dt & (42) 

Another representation for matrices with non-zero trace will be given in 
the next section (see eq. (|B5D). 



Finally, let us derive a kind of representation (40) for the case 



M 



B{t) = Y,C 3 {t)®D 3 {t) 
where Cj and D 3 are N\ x JVj and N 2 x ./V2 matrices. Let a 



(43) 



1,2 Ti = 1, 2, . . . , N{ be two sets of annihilation and creation operators that 
act in some Fock space, and n 101 is a projector on the space 



Hi 



{^1)^(2) J°>} 



(44) 



One can write 



; T-J i / dt B 

tr r e J o 

= trlli^i P exp ji^ rft^(a} 1) C i a(i ) )(a[ 2) D i a (2 ))| (45) 



The last formula ia an analog of (p~9|) . Then, repeating the calculations that 
have been done for derivation of eq. (^) from eq. (IB]), one obtains: 



2 ' 



tiPe'fi* 3 = [D 2 z 1 D 2 z 2 l[5(zl i) z ii) -l 
exp jz / ^4)^) 

I Jo L i=\ 
M 1 

+ E( Z (i) C W) - 9 trC i)(42) D i a (2) 



Zra - - tr D„ 



(46) 



Similar arguments allow to obtain the analog of eqs. (|38|) , (|40|) , and ( p]) 
for an matrix i? defined by eq. (f43|). 
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3 Alternative proof of bosonic path integral 
representation for the trace of path ordered 
exponent 

At first, we will evaluate the integral (^) assuming that B G su(N) . One 
can write: 

Z = lim / D 2 z5{( 1 dtJz-l--)e- e ti dtz " z e l ti dt ^ z+z " Bz) 
^+oJpbc Jo 2 J 



DC 



lim / ^ e -Mi+f) / D 2 ze jl d t z ^-± + iB + ^) z (4?) 

J 27T JPBC 



— oo 



Performing functional integration, one gets: 

Z = lim / P- 7 , (48) 

2tt det (-f t +i(B + V )-e) pBo 

To evaluate integral in (f48"D, let us consider eigenvalue problem 

^-A + z(£ + r7)-eW) = A<K*) (49) 

0(0) = 0(1) (50) 
The general solution of eq. fl4T)|) is 

= e^+^U^x (51) 

where Z7(i) is a solution of flT5j ) and vector x doesn't depend on t. 

Let e lOLr and £ r , r = l,...,iV be eigenvalues and eigenvectors of the 
matrix £7(1): 

U{l)i r = e ia ^ r (52) 
One notes that £7(1) e SU(N) and so a r are real and 

N 



J2a r = (53) 



r=l 
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Further, to satisfy ( |50"D we must put 



A = X r 



-e + ir\ + ia r + 2nin, n = 0, ±1, 



(54) 



and x — £r- So 



det -_ + i(B + „)- £ 



N oo 

J][ Y[ (27rin + ia r -Mr? — e) 



PBC r=l re=-oo 



oo 

-i ]~J (27rin) 



n=— oo J r=l 



r jv 



H(a r + ri + ie) ]} 1 + 



oo 

-i JJ (27rin) 

n=— oo 



r TV 



Y[(a r + ri + ie) J] 1 



« r + 7/ + ze 
27m 

(a r + 77 + ie) 
47r 2 n 2 



(55) 



r=l n=l 

Then, omitting irrelevant infinite constant and using well-known formula 

,2\ ! 



n(i-^]=£ 



sin rca 



n=l 



7r / 7ra 



one gets 



det (--^ + i(B + 77) - e 



AT 

n 



sin 



aj r + 77 + ie 



/ P_BC r=l 

Substituting ( |57|) in (f4~8"D, one obtains contour integral 



(56) 



(57) 



-^1+4 



Z= Hm / 7 



(58) 



One can close the contour of integration in Q58D in the lower half plane 
and represent Z as the sum of residues in the poles 



Vrn = —Ctr ~ i€ + 2tTTI 

The contribution of the pole with given r and n is 

e i(a r +ie) 



2ttz(-1) 



n+l ^ 

I] sin(^f^) 

l<s<N ^ ' 



(59) 



(60) 



s^r 
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So, omitting again inessential constant, one gets: 



lim 

r,n \ 



contribution of the residue in 
Vm = —&r — ie + 2nn 



= Et^" (61) 

r=1 II sin (^y^) 

Finally, using elementary but rather non-obvious identity 

Elv^- - = (-2i) lf - 1 \e^ £e- (62) 



' N / \ 
i rj sin 



r=l 



(that is valid for any complex numbers Appendix A for proof), one 

obtains: 

N j 

Z = J2 e iar = tr C/( 1) = tr P e* Jo dt m (63) 

r=l 



This proves the representation (|39|). 

Throughout the proof we didn't control inessential numerical normaliza- 
tion factors arising in front of functional integrals, determinants, etc. They 
can be easily reconstructed from normalization condition 

Z = N (64) 

B=0 

The representation (|39| ) is valid for traceless matrix B. In general case 
the following representation is valid 



Z = e-U> to ^) / D 2 z 6 (Y dtzh - 1 - f ) e'/o^Wj,,, (g5) 



This representation is alternative to one given by eq. (|42[). 

To prove (|65|), it is sufficient to repeat the proof given for representation 
( |39|) but without using the condition (|53l) . One can trace the cancellation 
of the pre-integral factor in (|65|) and contribution of the factor 
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JV 

e r=1 

Now let us check the representation fl40|). Representing 5-function in (fi"0|) 

as 

l[5 (Jz -1-^=Jd v e'JX* 1 - 1 -*)" (66) 
and integrating over 2^, z, one obtains the analog of (f48[): 

Z = lim / £ty(t) y— r (67) 

^+°J det(-| + i(B(t)+T)(t))-e} 

The eigenvalue problem 



PBC 



'-j t + + - = A0(t) (68) 

0(0) = 0(1) (69) 



has the solution 



A = \ rn = — e + ia r + i dtrj(t) + 2mn, n = 0, ±1, . . . (70) 

Jo 

where numbers a r are defined by equation Q5^). So determinant in flB7|) 
depends on ^(t) only via 

V=[ 1 dt V (t) (71) 
Therefore, if one introduce new variables 77 and 

77 n = f'dte^vit), njtO (72) 
Jo 



instead of in (|67|) , one finds that integration with respect to r\ n gives only 
inessential constant and so the integral fl6"7|) transforms just in the integral 
([48]). But this means that we reduced the proof of representation ( |40"1) to one 



of the representation (1391) which we have already proved. 
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Finally, let us prove the validity of the representations (|38| ) and fl4T|). The 
proof can be reduced again to one of the representation ( |39D by means of the 
formula 



lim det (-4- + iBe~^) = e*/o dttrB ^det (~ + iB(t)) 

\ dt J \ dt J 



(73) 



that is valid for any N x N matrix B(t) (the proof of ( |73[ ) is given in Ap- 
pendix B). Indeed, evaluating the integral with respect to z, z in (|38|) , one 
obtains, after inserting of the factor exp (-6fodtz^(t)z(t + e) 



dfj 6 
lim lim / -, j-r (74) 



5 ^ + o e ->+o 7 2vr det (-| + i(S + 77 + z5) e - £ i 



-co \ a!t 

and further application of (^) reduce (|74D to (|48|) . (Remind, that tri? = 
by assumption.) 

The same arguments allow to check the represent at ion (|4T|) . 

Thus we have checked all four representations for the trace of path ordered 
exponent derived in the section 2. We see the proofs given in the present 
section are very unlike ones given in the section 2. This can be considered 
as strong confirmation of the validity of the representations obtained. 



4 A variant of non-Abelian Stokes theorem 

Analogs of Stokes theorem for Wilson loop 

Z = ^Yi^ dx " A ^ x) (75) 



were proposed by several authors |^3) . In this works Wilson loop is expressed 



via some area integral over surface spanned on 7 with rather complicated path 
ordered prescriptions. Recently Dyakonov and Petrov derived another, very 
smart variant of non-Abelian Stokes theorem. Their result we already cited 
(see©). 

In this section we will prove a new variant of non-Abelian Stokes theorem 
for G su(N) that is similar to one given by Dyakonov and Petrov. For the 
case iV = 2 we will be able to derive from our results the Dyakonov-Petrov's 
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formula (|H]) but in slightly corrected form. The little discrepancy between 



our results and those by Dyakonov and Petrov arises, most likely, because of 
some subtleties in the definition of the corresponding functional integrals. 
Let the path 7 in (|75|) be parametrized as 

7 = { x » = q ^t), < t < 1, g"(0) = <?"(*)} (76) 

Then 

Z = tlFe ifodtr(t)A^q(t)) ^ 

and we can apply the representation (|40|): 

Z = J D 2 ^z - 1 - ^) e fo dtzt (-£^ A ) z (78) 
Let £ r = £ r (x), r = 1, . . . , N be any field such that 

T(g(t)) = S{t) (79) 

The formula ( [78"D can be rewritten as 

Z = Jd 2 £ It 8 (?(x)t(x) - 1 - f ) e-^ tD ^ (80) 

where .D^ = <9 M — iA^ is the usual covariant derivative. Further, applying the 
classical Stokes theorem, one can easy to prove that 

j> dx» gD&x) = / s dx» A dx v (d^D^ - l -S}F^ (81) 

where E is any surface for which <9E = 7. 

Substituting (|8l|) in (|80|) , we get our variant of non-Abelian Stokes theo- 
rem: 

Z= trPe'^"^ 

N 



(82) 



Now let us transform eq. (^) into the form that would be similar to eq. 

(ED- 
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Let A a , a = 1,...,N 2 — 1 be Hermitian generators of SU(N) in the 
fundamental representation. They can be normed as 

trA a A 6 = 25 ab (83) 

and satisfy equations 

[A a , \ b ] = 2if abc \ c (84) 

[A a , \ b } + = ^S ab I N + 2d abc X c (85) 

One notes that 

A, = y a I F ^ = = y (<V a , ~ d »K + f ahc K A D (86) 

Matrices A a also obey "Fierz" identities: 



ki + N 

if abC ^ij^kl^mn = 2(fiilfimjfikrn ~ ^in^jk^ml) (88) 



2 KjKi + j^hjdki = Su5 kj (87) 



One introduces new variables 



I a (x) = ?(x)\ a Z(x) (89) 

Variables I a are not independent. Using "Fierz" identity (|87j), formula 
9|) can be represented as 



Let / = I a X a , £t£ = c. One notes that c = 1 + y on the surface of 
integration in fl82|) . 

The matrix / can be represented in the form 

I = Udiag(ai, . . . , cln)U^ (91) 

where 
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$>i = 0, UeSU(N) (92) 
i=i 

The matrix in R.H.S. of eq. (j90|) is a projector on the one dimensional 
subspace. This means that eigenvalues of the matrix in the L.H.S. of eq. (|90|) 
are all equal to zero except only one that is equal to 1. So, up to renumbering 
of eigenvalues, 



2c 

Oi = -J7, Z=1,...,JV-1 

a N = 2c (iV ~ X) (93) 



One notes that 



2c 



JV 



(2c) fe (JV- 1) 
JV* 



(JV 



\k-l 



-If 



(94) 



Thus the matrix / has — 1 coincident eigenvalues. But this means that 
the solutions of eq. (0) are in one-to-one correspondence with the points of 
a coset SU{N)/U{1) x SU(N - l)w CP N ~\ 

Indeed, the matrix U in formula (^1|) can be presented in the form 



U = U 1 V 1 V 2 



(95) 



where 



Vi 







, V x G SU(N - 1) 



(96) 



V o-.-o i J 

V 2 = exp{diag(t, ...,t, -(JV - l)t} 



and a matrix C/j represents an element of the coset SU(N)/U(1) x SU(N-l). 
So, due to coincidence of the first JV — 1 eigenvalues of the matrix /, one finds 
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/ 2c 

I = U\diag I - — , 



_2c 2c (iV-2)\ ^ 



(97) 



Therefore the set of solutions of eq. 
SU(N-l). 

Further, using identities fl3T|), 

c = const 



is isomorfic to coset SU(N)/U(1) x 
8f), one finds that on the surface 



where 



1 



2 tvID^ID^I 



- tr 

2 ^ 



D M J = 9 M /-i[A M ,il 
Now let us insert in (|82| ) an identity 



1 = n / n^w^o^o*) - /a ( a; )) 



By virtue of (051), (|9~9T), the eq. can be rewritten in the form 



(98) 
(99) 

(100) 
(101) 



Z = J Dfi(I) exp < -z J dx^ A dx 1 " 



tr ID\ U ID U ]I — - tr//, 

8c 2 [M 1 4 



/if 



where 



(102) 



D^i) = (u Dia ) J £>we - 1 - f ) - n 



(103) 



The consideration given above shows that Dfi(I) is nothing but SU(N)- 
invariant measure on the coset SU(N)/U(1) x SU(N — 1) 

The explicit form of the measure Dfi(I) is rather cumbersome but, for- 
tunately, it appears that measure Dfi(I) in ( |102j ) can be replaced by the 
measure 
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N 

Dfi'(I) = (DP) J] 6(ti I k - c KN ) (104) 

k=2 

where numbers Ck,N are defined by formula (0) in which one must put 2c = 
iV + 2 by virtue of the first 5-function in eq. (|103|) . 

Indeed, numbers tr I h define uniquely characteristic equation for the ma- 
trix / and, consequently, its eigenvalues. This leads to representation (|97|). 
So on the surface tr I k = c^^n any functional is invariant under trans- 
formation 

$(/) -> $(ViV£IVjV?) (105) 
where a matrix V1V2 is defined by eqs. (|96|). Therefore in any integral 

J DfiVWI) (106) 

one can perform all integrations except those corresponding to integration 
over coset SU(N)/U(1) x SU(N - 1): 



J Dfj?(T)$(I) = const J Dn(I)$(Uidiag(-^, 



iV' 2c N )Ul) 



(107) 

(see (|97D). But this just means, in particular, that one can replace the 
measure Dfi(I) in ( |102| ) by the measure Da' (I) defined by eq. ( |104| ). 
Now we can formulate our final result: 



tr P exp li j> dx^ A 

N 

= DIl[5(tTl k -c k 

J k=2 

exp < i J dx^ A dx v 



N 



1 



2(iV + 2 



-tr(ID^IDJ) + ^tr/F Mi/ 



(108) 



where 
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^-T^V -^' + t-l)'] (X09) 

Let us compare our results with those due to Dyakonov and Petrov [gj]] . 
Putting in (fTUS]) N = 2, I = I a a a and the changing J a -» 2/ a , we get: 



tr P exp jz ^ dx^ 
= J DI5(I 2 - 1) exp |~ y dx" A rfx" (e abc I a D^ b DJ c - / a F;,) j 



(110) 



Comparing formulae (|110|) and fll4f ), we see that they differ by the fac- 



tor 1/2 in front of the "action". Most likely this discrepancy arises because 
of some subtleties in the definition of the functional integrals that play the 
important role in our discussion. In particular, if one ignores the difference 
between functional integrals constructed by means of Weyl and normal sym- 
bols, one obtains the following representation instead of 



Z = f d\ n 5 (e f (x)e(x) - 1) e -§^$ D *w (in) 



Further, using formulae (|98| ) and (|9^ ) with c = 1 instead of c = 1 + ^> one 
can easy trace that representation (|1 1 1|) leads exactly to Dyakonov-Petrov 



formula (0). But representation ( |111| ) is wrong. So just formula ( |110| ) must 



be considered as correct version of non-Abelian Stokes theorem for N = 2. 

This doesn't mean, however, that Dyakonov-Petrov formula ([Tip is incor- 
rect. But it means that the definition of functional integral in (|T4] ) must be 
clarified. 

Another discrepancy between our results and those due to Dyakonov and 
Petrov arises in the case N > 3. Dyakonov and Petrov pointed out in 
their work [^T| that integration in the functional integral representation for 
Wilson loop must be performed over the coset SU(N)/ [U(N)] N_1 whereas 
in our representation ( |108| ) integration is carried out, in fact, over the coset 
SU(N) /U{1) xSU(N— 1). But nowadays it is hard to discuss this discrepancy 
because no explicit formulae for the case N > 3 were given in ||21||. 
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5 Bosonic worldline path integral represen- 
tation for fermionic determinants and Green 
functions in Euclidean space 

Bosonic worldline path integral representation for fermionic determi- 
nants can be obtained directly from formulae fllCf ) and results of the 
section 2: 



In det(zV + im 
1 r°° dT 



2 I" T ^ La D * DW " 1 - T) 4 - 3 ) 
At exp | J dt 

-T(^a^)(z^F^z 



q 



(112) 



Here variables z = {z r , r = 1, . . . , N} and ip = i = 1, , 2, 3, 4} 
describe colour and spin degrees of freedom respectively, z' and i[)> are 
complex conjugated to z and ■0, 



N 



° 2z = UY[d(Rez r (t))d(Imz r (t)) 

t r=l 
N 

D2 ^ = II II d(Re#(*)) rf(Im#'(*)) 



(113) 
(114) 



t r=l 



and At is a normalization constant. The latter can be evaluated from the 
condition 



<x|tre- T (- v " vw ^b> 



= 4N < x\e T9 ^\x >= ^ - 
L4=o ' ' 2tt 2 T 2 



(115) 



Indeed, putting A = in (|1121) and comparing the result with ( |115|) 
one obtains: 
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jVn 



-1 



2tt 2 T 2 



A X ' / Dq [ D 2 i,D 2 z5(zh-l-- 

J JPBC \ 2 / 

<?(0)=<?(l)=z 



5 V - 3) exp I J dt 



q 2 

AT y v 



(116) 



Obviously, Mr doesn't depend on x. 

Remind, that in Euclidean space ib'ib is SO (A) scalar. So representation 
( |112| ) is manifestly relativistic and gauge invariant. 

Our next task is to derive bosonic worldline path integral representa- 
tion for Euclidean Green functions. In what follows, we restrict ourselves 
only to derivation of such representation for generating functional Z(j) for 
vacuum correlators 

< ip} i (x 1 )iP fl (x 1 )...ib} n (x n )iP fn (x n ) > (117) 

However, our method is quite general and can be easily applied to deriva- 
tion of analogues representations for arbitrary Green functions. 

Standard functional integral representation for Z(j) can be written as 

ZU) f I)A ( ^^ J P^Umi^- 1 ' 11 - ■"'^-''' ! ^'^ (118) 

Integrating with respect to fermionic fields, we get 

Z(j) = f DA e SvM J] det(zV + im f - ij f ) (119) 
/ 

So our task is reduced to derivation of bosonic path integral represen- 
tation for determinant 

det(iV + im — ij) 

The latter problem can be easily solved by applying of the results 
obtained in the section 2. Indeed, 
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lndet(zV + im — ij) = - ln[det(zV + im — ij)^ 5 ] 2 
= ilndet(V M V -a^F^ + dj + (m-jf) 

2 Jo T 

Using representation for path ordered exponent from section 2, one 
obtains: 



In det(zV + i(m + j)) 

i e _ m2T r DqD 2^ D 2 z5 (j z _ 1 _n\ s 

Jo T Jpbc V 2 / v 



2 



<? 2 



exp j y (it 

-T^^)(z i F^z-^^d^q)-2mj(q) +f(q)]} (121) 



Substituting ( |120| ) in ( |121| ) for each /, we get worldline bosonic path 



integral representation for generating functional Z(j). The corresponding 
formulae for n-point correlators are rather cumbersome but quite com- 
putable. Author hopes that they can be used for computer simulations on 
the lattice. 

6 Bosonic worldline path integral represen- 
tation for fermionic determinants and Green 
functions in Minkowski space and quasi- 
classical approximation in QCD 

The derivation of bosonic worldline path integral representation for 
fermionic determinants in Minkowski space is slightly more involved then 
one in Euclidean space. The origin of complications is non-unitarity of 
finite dimensional representation of the Lorentz group. 

The analog of the representation @ in Minkowski space can be written 

as 
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lndet(*V - m) = -- [°° % e ~ tm2T tr e iT ^ v ^ + ^ F ^ (122) 
2 Jo T 



Trace in (122) can be represented as functional integral: 



nv^-v.v^-^/.;,.,, = / Dq tr p^/o^^^HT^w) (123) 



PBC 

eq. ( |123| ) is an analog of eq. flHiD. Path ordering in ( |123| ) corresponds to 
colour and spinor structures. 

However, in contrast to Euclidean case, we cannot directly use the 
representations of the type (fH]); ( f46[ ) to write ordered exponent in ( |123| ) 
as functional integral. Indeed, those representation comprise, in particular, 
the factor 

5(^ip - 3) (124) 

(see ( |112p ) that is not Lorentz invariant because spinor representations of 
Lorentz group are not unitary. 

To obtain manifestly Lorentz invariant representation, we will use, at 
first, representation ([4(]) for describing of colour degrees of freedom and 



representation (41) for describing of spinor ones. In such terms eq. ([12 



can be rewritten as 



tre i(-V M W+(r'«'j} Bf ) 

r N 
= lim / DqD 2 zD^D^DX5(z^z - 1 ) 

J PBC 2 



JVt exp 



dt 



q 2 

h iz'z + iibib+ 

AT ^ 



z ] qA{q)z + T^a^e^i^z^F^z) + A(^e" e % - 1)] } (125) 



In eq. ( 125 ) the measure D 2 z is defined by ( 113 ) but ip and t/> are 
independent complex variables. Mr is a normalization constant that will 
be computed later. 

Integrating over if), if) in (|125|), one gets: 
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r N 
= lim / DqD 2 zD\5(z ] z - 1 * 



JVrdet -1 
exp < % J dt 



+ iTa^(z^F uv z)e~ e tt + iAe" e * 



(f 



+ i-z^i + z*qA(q)z - A 



But 



(126) 



lim det 1 



d 
It 



= e- 2l ti dtm det- 
by virtue of identity ([731). Further 



-^L+iTa^^F^z) /A 
at 



(127) 



det" 1 



dt 



det" 1 



- 7 ^ + 7 zT^V^) + *7°A 



128) 



The operator in R.H.S. of ( |128| ) is anti-Hermitean. So we can write 



det" 1 



+ -fiT^^F^z) + 2 7 °A 



P_BC 



D 2 1p€ L fo ^(# t 7 V'+T(V t 7 ^^)(^ t i^^+A^WV) ^29) 



In the last formula ^ and ip are already complex conjugate variables, 
measure D 2 ip is defined by eq. ( |H4| ), and the "action" is real. So functional 
integral ( |129| ) is well defined. 

Introducing standard notations tp = tp^ and substituting ( |127| )- (|129|) 
in (|126|), one obtains, after integration with respect to A, 
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= J DqD 2 i)D 2 z5 (z*z - 1 - y ) 5 
Mr exp {if dt 



- 3 



z^qA{q)z + T(i/ja^i))(z^F^(q)z)\ } 



(130) 



The normalization constant Mr can be computed in the same way as 
its analog in eq. ( |112j ): 



M T l 



N 



2ir 2 T 2 



I Dq [ D 2 ^D 2 z5 (zh - 1 

J JPBC V 



q(0)=q(l)=x 



N 

~2~ 



— 3 J exp < i J dt 



q 

— — + iz'z + iip^ip 



(131) 



Substituting (|130|) in ( 122 ), we obtain desired representation for fermi- 
onic determinant: 



In det(zV — m) 
1 r°° dT 



-im^T 



2 Jo T 
Mr exp < 



N 



dt 



DqD 2 tpD 2 z5 [zh-1- 
pbc V 2 

q 2 

h iz'z + ixbib+ 

AT vv 



5 - 3 



z*qA(q)z + T(i)a^)(z^F fJ ' u (q)z)] } 



(132) 



The representation ( |132| ) is manifestly gauge and Lorentz invariant and 
comprises only bosonic variables. The "action" in ( |132| ) is real. So, as we 
will see soon, it is convenient for application of stationary phase method. 

Now let us derive bosonic path integral representation for generating 
functional Z(j) of gauge invariant Green functions 



Gf 1 ,..., fn (xi,...,x n ) =<T(i) h (x 1 )i) h (x l )...i) fn (x n )i) fn (x n )) > (133) 
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The derivation is completely analogous to one given in the previous 
section for corresponding Euclidean correlators. 

For Z(j) there exist standard path integral representation via anti- 
commuting variables: 



Z(j) = / DAe / D^D^e f 

= f DAe iSrM J] det(iV -m f + j f ) (134) 
J f 

Repeating with minor changes the derivation of ( |121| ), one gets 



In det(iV — m — j) 

- -\ r % e " m2T La d « dw - 1 - t ) s - 3 

Mt exp < i I dt h iz' ' z + iipip + z^qA(q)z+ 

[Jo [AT 

T(^^F^(q)z) + 2mTj(q) - iT^fW^q) - Tf(q)] } 



(135) 
2 

Substituting ( [L35f) in ( |134j) , we obtain worldline pat integral represen- 
tation for Z(j). 

Our next task is the investigation of quasiclassical approximation in 
QCD. To this end, we will formulate a scheme of evaluation of two-point 
function 

Gf (x,y) =< T(^ /o (x)^/ (x)^/ (?/)^/o(l/)) > (136) 

This scheme can be easily generalized for evaluation of arbitrary Green 
functions. 

The equations for the stationary point will be interpreted as quasiclas- 
sical equations in QCD. They will be formulated in terms of particles that 
have spin and colour degrees of freedom and interecting with Yang-Mills 
field. 
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First of all, we introduce more condenced notations: 



Qf = {9/,4' z /'^/'^/' T /' m /} ( 137 ) 

D Qfi ...) E -iff J pBc n qDW ^ z -i-^ 



S[Q f ,A] 



2 jo 

5 - 3) 7V T (- • • 
1 



eft 







AT 



+T f (iJ f a^ f )(z}F flv (q)z f ) - T f m) 
Further, using eq. ( |135| ), one can get: 



(138) 



(139) 



In det(iV — rrif + j 



j=0 



DQ f e iSlQ f' A] R(x\Q f ) (140) 



S 2 



In det(iV — rrif + j 



j=0 



Sj(x)5j(y) 
= J DQ f e iS ^R(x\Q f )R(y\Q f ) 

+ i5{x-y) J DQ f e iS[Q f' A] T f £ dt x dt 2 5{q f {h) 



-Qffo)) (141) 



where 



R(x\Q f ) = T / dt[2m f - iMQ-fMWM* ~ ?/(*)) (142) 
For any functional W(j) 



W(j) 



5 2 \nW{j) 5lnW(j)ShiW(j) 
Sj(x)8j(y) 5j(x) 6j(y) 



Using (pi, O, (HU) and applying O for Wtf) 
m/ + j), one obtains: 



(143) 
det(iV - 
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G f0 (x -y) = GfJ(x -y) + Gf>(x - y) (144) 

where G$ and G^ correspond to the first and to the second terms in 
R.H.S. of O respectively: 

G^(x-y) = J DADQ f0 R(x,Q f0 )R(y,Q f0 ) 

exp liS YM + iS[Q fo , Aj+^jDQj e^ Q f^ 

+5(x-y)(---) (145) 
Gf{x-y) = J DADQ fo DQ' fQ [R(x,Q fo )R(y,Q fo ) 

exp jzSVw + iS[Q fo ,A] + iS[Q' fo ,A] 

+ E/ DQ f e tS[Q f> A] ^ (146) 

In (HD (• ■ ■) means the factor at S(x - y) in R.H.S. of (fl41]) . 
At first, we investigate the function G^. 

The function Gf Q (x — y), in itself, is defined up to counterterm 

const 5 (x — y) 

by virtue of ultraviolet divergences. Then the last term in R.H.S. of (|145|) 
only redefines this counterterm and so can be omitted. 
Expanding 

exp jDQ f e lS ^ (147) 
in series, we can represent (|145|) as 
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Eff^TT / DADQfy (n fl DQif) R(x,Q f0 )R(y,Q f0 ) 
n f 11/™/! •> y f jf=1 J 

{n f 
iSym +iS[Q fo , A] + E iS[Qj f ,A] 
f i/=i 

E^(^!/;W}) (148) 



7o 

n f 



Obviously, each function Gf^{x — y; {nf}) correspond to contribution 
of all diagrams comprising nj 1 quark loops of flavor fx, rif 2 quark loops of 
flavor / 2 , etc. 

We will investigate each term in the series ( 148|) separately. At first we 
consider the term with nf = 0: 



G^ix - y; {n f = 0}) = J DADQ h e iS ™ +iS ®f°^R(x,Q fo )R(y,Q fo ) 



(149) 



/o 

The function R(x\Qf) can be represented as 

R(x\Q f ) = J d 4 p J 1 dtR{p,t\Q f )e ip{x ~ q ^ )) (150) 

where 

R(p,t\Q f )=T(2m f + p^(t)Ym) (151) 
Then we change variables: 

q(t)^q'(t) = q(t)-q , A{x) -> A{x - q ) (152) 
where the function q'{t) obeys the boundary conditions 

q'(0) = q'{l) = (153) 
By virtue of translational invariance we get 

= I DA I DQ f „ I dq e iSYM+iS[Q f0 ,A] 

J Jq(0)=q(l)=0 J 
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1 rl 



d pid p 2 dti dt 2 R(p 1 ,t 1 \q fo )R(p 2 ,t 2 \qf ) 

e ipi(y-qf (ti)) e ip2(y-qf (t 2 )) e -i(pi+P2)qo 

(154) 

Integration over q gives 8{p% +p 2 ), and we obtain the following repre- 
sentation for Fourier transformation of G^(x — y; {rif = 0}): 



dxe- ipx G ( l\x;{n f = 0}) 

1 dh f 1 dt 2 f DA f DQ h R{pM\Qh)R{-pMQh) 

JO J Jq fo (0)=q fo (l)=0 



exp{iS YM + iS[Q fo ,A] - ip(q fo (ti) - q fo {t 2 ))} 



(155) 



Now it is already easy to write equations for stationary point for the 
action in ( |155| ). Omitting the index f , we get: 



\v u F a ^ = f 1 dtl a (t)q"(t)5(x-q(t)) 



dtP^S^Six-qit)) 



i (j t - iq*Ar(q)\ z + TS^F^z = 



dt M 



-q, + <fF;j a + V„J? p ( ff ) W = p^fa) - 5(t 2 )) 



l r l 



"p 

l r 1 



AT 2 



( dtq 2 + - [ dtS» u F a T a = m 2 
Jo 2 Jo * 



(156) 

(157) 
(158) 
(159) 
(160) 



where 



I a = z^X a z, = tPct^tP 



(161) 
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Eqs. ( |156| )-( |I~59|) can be derived by variation of the action in ( |155| ) with 
respect to A®, z, ip, and q p . In the derivation of ( |159| ) we used (|156|) - (|158| ). 
The eq. ( |16U| ) is obtained by differentiation with respect to T. 

It easy to obtain closed system of equations in terms of A^, z, I a , and 



Let 



Then 



where 



\ a I a (162) 



iDI = TS^[F^,I] (163) 

dgp, = 2TPF a[fl n S u]p (164) 
dt p y ' 



DI =^_ i[rMq)J] (165) 



dt 

Unknown functions in (|156|)-(|ToTS|) also obey boundary conditions 



g(0)=g(l) = 0, z(0) = *(l), i/,(0) = ^(1) (166) 
They also satisfy the equations 

N 

z x z = 1 + —, ijjijj = 3 (167) 

because of the presence of ^-functions in the definition ( |138| ) of the mea- 
sure DQ. Apropos, we didn't introduce Lagrange multipliers to take into 
account these (^-functions because the solutions of eqs. ( |157| ), ( |158[ ) auto- 
matically satisfy the conditions 

z* z = const, ijjtjj = const (168) 

Equations Q156Q , ( |159| ), and ( \L63j) are nothing but generalization of 
well-known Wong's equations p4| that describe classical spinless particle 
interacting with Yang-Mills field. Indeed, if one omits the terms containing 
the tensor of spin S pu in eqs. (|156|) , ( |159| ), and (p. 63) one gets just Wong's 
equations up to term 
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-p^O-W) (169) 

Eqs. (|156|) -( |I"58| ) admit simple interpretation. At "time" t\ quark- 
anti-quark pair with momentum p is created. Then quark and anti-quark 
move interacting with Yang-Mills field. The union of quark and anti-quark 
trajectories forms a closed loop passing the point q = 0. (See (|166| )). At 
the "time" ti the quark-anti-quark pair annihilates. 

An analog of eqs. ( |156D -( |ToDD for terms with nj ^ in expansion (148) 



can be derived in the similar way. Let S^ nf=0 ^ be the action in the formula 
(|155|) and J afl [Q} be the current in R.H.S. of eq. (|156|) . Then the analog 
of S^ n f =0 ^ for the term with non zero numbers {nf} in ( 148 ) is 



S {n f } = S {n f =0} S[Q 3r A] (170) 

/ j/=i 

So instead of (|156|) we have in this case equations 



iv„iw = J a nQ f0 ] + e E ja nQ jf ] (i7i) 
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Equations 



«™ 0, f^ = 0, (172) 
OZj, dl> 

coincide in form with (|158|) , ( p.60|) , and (|161|) , whereas the equation 

SS {n f } 

= (173) 

5q jf 

differs from (|159|) by the term (|169|) . 

We see that quasiclassical configurations that give the main contribu- 
tion in functional integrals (|148 ) are defined by equations of very similar 



structure. The same statement is valid for semiclassical configuration that 
give the main contribution in the function G^ 2 '(x — y) as well as in any 
other gauge invariant Green functions. Therefore proposed scheme seems 
to be sufficiently general for application in QCD and other gauge theories. 

In this paper we restrict ourselves to formulation of general scheme 
of quasiclassical approximation in QCD leaving elaborating of details as 
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well as applications for forthcoming papers. So at this point we stop our 
investigation of quasiclassical approximation in QCD. Brief discussion of 
possible applications will be given in the next section. 

7 Conclusion 

In the present paper we derived, first, new path integral representations 
for path ordered exponent (see eqs. (]38D-(]4lD, fl46|) , (|65|)). We give 

two alternative, entirely independent derivation of these representations. 
So these results seems quite reliable. 

Then we applied these representations to derive new variant of non- 
Abelian Stokes theorem. Our result is represented by formula ([82]). Then 
we transformed the latter to obtain another formulation of non-Abelian 
Stokes theorem that is similar to one proposed recently by Dyakonov and 
Petrov I^TJ. As a result, we got corrected and generalized version of the 



theorem proved in [21 



Dyakonov-Petrov version of non-Abelian Stokes theorem was already 
used in discussion of the role of monopoles in QCD |21| and in attempts 
to derive string-like effective action in framework of QCD |^5 |. So one can 



hope that our more general and simple version of this theorem will be also 
useful in discussion of various problems of QCD. 

Further, we derived pure bosonic worldline path integral representa- 
tions for fermionic determinants as well as fermionic Green functions in 
Euclidean QCD. (See eqs. (|T|), (|120D , and fll2lD ). 

This representations comprise only integrations with respect to bosonic 
variables. On a finite lattice all integrals are quite simple and well conver- 
gent. (Remind, that domain of integration with respect to z and ip in eqs. 
( |112| ), ( |120|) , and ( |121| ) is bounded). So representation derived seem quite 
appropriate for lattice simulations. 

Our results for fermionic determinant and Green functions can be used 
also in another way. 

Namely, if one substitute instead of Wilson loop some phenomenologi- 
cal anzatz, one obtain formulation of the theory purely in terms of point 
particles. The most well-known example of such anzatz is Wilson area law: 



< tr P exp j> dx^ >= exp < —KS min + 



perturbative 
corrections 



(174) 
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This anzatz was applied, in particular, to derivation from QCD a 
quark-anti- quark potential used in potential models (see, for instance, 
recent papers PB[ , |27| and references therein). 

Other anzatz for Wilson loops were proposed in framework of Dual 
QCD model |28| and stochastic vacuum model [^]. (See also [^] for 
comparison of results obtained in framework of these models). Recently 
a sting-like expression in spirit of stochastic vacuum model was derived 
in papers [p5|| , |30]. Anzatz of another type, that gives an expression for 



Wilson loop in terms of trajectories of monopoles, was proposed in [^T 



All these results could be combined with ours to derive some effective 
action in terms of point particles that correctly describe colour and spin 
properties of quarks. Obviously, that such theory is much more simpler 
than initial quantum field theory and so thus approach of investigation 
seems to be rather perspective. 

Finally, in the present paper we also derived bosonic worldline path 
integral representation for fermionic determinants and Green functions in 
Minkowski space and started the investigation of the quasiclassical approx- 
imation in QCD. 

The key point in the formulation of quasiclassical approximation is 
quasiclassical QCD equations (|156|) -( |l60l) and (|170|) -( |173|) which arise nat- 
urally when one applies the stationary point method to evaluation of func- 
tional integrals that defines Green functions in QCD. Quasiclassical equa- 
tions derived appear to be nothing but a generalization of well-known 
Wong's equations. 

We formulated only those quasiclassical equations which arise in the 



problem of evaluation of the concrete Green function ( |136|) . However, our 
method is quite general and one can easy to derive analogous equations in 
generic case. 

The next problem in investigation of quasiclassical approximation in 
QCD is the solution of quasiclassical equation of motion. Though these 
equations are very complicated, this problem doesn't seem hopeless, at 
least, in non-relativistic approximation. In electrodynamics in non-relati- 
vistic approximation one can neglect bremsstrahlung and retarded effects 
and, as a result, reformulate initial theory in terms of particles interacting 
by means of Coulomb forces. In the same way one may hope to derive po- 
tential of interaction of heavy quarks in QCD. Such approach is alternative 
to ones based on various anzatz for Wilson loops. 
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Another interesting possibility to understand the quark confinement 
in framework of quasiclassical approach is connected with existence of 
classical solutions of Yang-Mills equations with singularity on the sphere. 
Such solution were discussed in the context of the problem of confinement 

though the existence of such solutions was 
Some solutions with singularity 



recently in the papers [0T 



pointed out by several authors in 70's [32 



on the torus and cylinder was discussed in papers |33 



If solutions with singularity on closed spacelike surface existed also for 
eqs. ( 156|) - (|160 ) and their modifications, then quarks, moving inside such 
surface and interacting with Yang-Mills field, couldn't cross the surface. 
This would mean the confinement of quarks. So it is interesting to in- 
vestigate singular solutions of equations ( |156D -( |ToT]| ) and, if such solutions 
exist, to develop the corresponding quasiclassical perturbative theory. 



Appendix A 

In this appendix we shall prove the validity of the identity d62|) . 
We will prove fl62|) by induction. For iV = 2 it is easy to prove fl62|) by 
direct calculations because both sides of (^) are simple rational functions 
of variables exp{|a r }. 
Let us denote 

x r = e iar (Al) 
Then L.H.S. of fl62|) can be represented as 

JV 

II y/x r n 

i-^f- 1 n r=1 rr _ -r E(-l) i+1 *f n i*r-x.) (A2) 

l<p<q<N ' l<r<s<N 

whereas R.H.S. of (p2]) as 

N \ N 

fe=l / 7 = 1 



'-2i) N - 1 



\k=i J j 

So eq. fl62D is equivalent to algebraic identity 
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N I N 



j:(-iy +i xf n (sr-*,)= E** n (^-^) (^ 4 ) 

j=l r, B & \j=l I l< r <s<N 

l<r<s<N x / — — 

To prove (A4), we check, at first, that L.H.S. of (A4) is vanished if 
X{ = Xj. Obviously, it is sufficient to consider the case x\ = x 2 = x. 

If Xi = X2 then only j = 1 and j = 2 terms survive in L.H.S. of (A4). 
The j = 1 term is 



2<r<s<7V \p=3 / \3<p<g<iV 



whereas the j = 2 term is equal to 



II (xr-x s ) = -x N [ll(x-x p ))( n (^-^)) (^6) 

r,s#2 \p=3 / \l<p<g<JV / 

l<r<s<JV X ' N ~~ 7 

So terms with j = 1 and j = 2 are cancelled. 
Thus L.H.S. of (A4) can be represented as 

P(xi, . . . ,x N ) J[ (x r -x s ) (A7) 

l<r<s<N 

We must prove that 

N 

P(x 1 ,...,x N ) = '^2x r (A8) 

r=l 

L.H.S. of eq. (A4) is polynomial of degree N with respect to each 
variable Xj. Consequently, polynomial P(x±, . . . ,xn) is linear in each Xj. 
So 

P(x 1 , ...,x N ) = a(x 2 , x N )x 1 + b(x 2 , ...,x n ) (A9) 
Comparing coefficient at x^ in L.H.S. of (A4) and (A7), one find that 
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a(x, . . . , xn) = 1 (A10) 
and so it remains to prove that 

JV 

i=2 

To prove (All), let us compare the sums of Xi-independent terms in 
L.H.S. of (A4) and in (A7). They can be written as 

N 

E(-!) J+1 < II (-*.■) II (Xr-Xs) (A12) 



3 

j = 2 2<r<N 2<r<s<N 



and 



b(x 2 ,...,X N ) Y[ (- X r) II (Xr-Xs) (A13) 
2<r<N 2<r<s<N 



respectively. 

So (All) is equivalent to 



N N 



E(-i)^- 1 n {xr-x s ) = \Y J xA n ( 

j=2 r,s^j \j=2 I 2<r<s<N 

2<r<s<N V / - - 



But the latter equation is valid by virtue of induction assumption. 
Indeed, it can be obtained from (A4) by change N — > N — 1, , x\ — > 
x 2 , ■ ■ ■ , xn-i —* x N . This finishes the proof of identity d62|) . 

Appendix B 

In this appendix we prove the formula ([[3|). First, one notes that up to 
inessential factor 

del A* (fll) 

So it is sufficient to prove that 
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lim lndet \-^-e £ M + %B 
e^+o \ dt 



or, equivalently, 



_ kdet __ + iB 



PdttiB (52) 
j o 



.. . / d f ± 
lim trln — -e dt 
e-+o \ (it 



+ ijB ) ~ trln [^Jt +iB ) = ^Io dttlB 



Let ft is orthogonal projector in L 2 (S' 1 ) on any finite dimensional sub- 
space. Then 



lim 



trflm (-^e e ^ + iBj — tr ft In f-^ + i5 







(54) 



In particular, let 7ii be subspace in L 2 (S l ) that is orthogonal to one di- 
mensional subspace spanned on the function <f>o(t) = 1. Then 



lim 



trln [--^^ + ijB j _ trln + iB 



lim 



tr Wl In (-j/^ + iBj - tr Wl In (--^ + iS 



In the space Hi the operator 



d ,d_ 

— -e <u 
dt 



(55) 



(56) 



is invertible. Let G e (t — t') is the kernel of 



d d 

— e dt 
dt 



The function G e (t — t') can be expresseed in terms of eigenf unctions 

<P n (t) = e 2mn \ n^O (57) 
of the operator (B6) and its eigenvalues 
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\ n = -2irine 2nme , n ^ 



Obviously, 



_ (/> n (t)f n (t) 1 ~ sin 2im(t - f - e) 

Cr^t - t ) - 7 - -- 2^ 



n^O 

Using the formula 



7T 



n=l 



n 



(58) 



sin27rne 1 



71=1 



7m 



(that is valid for < e < 1), one finds that 



whereas 



So 



lim GJO) = - 



C o (0) = 



Go(0)^ £ lim o G e (0) 



But for -1 < i < 1, 0, 



e limG e (t)=Go(0 



(59) 



(510) 



(510') 



(511) 



(511') 



We will see soon that R.H.S. of (B2) is not vanished just by virtue of 
(BH). 

Further, up to inessential constant, 



tr Wl In (-^e^ + iB j = tr Wi hi 



1+i Hk*J B 



d 
"dt 



e dt 



B 



d d_ 

■—edt 
dt 



B 



d 
~dt 



+ 



(512) 
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By virtue of equation 

f 1 dtGJt-t') = (513) 
Jo 

one can replace tr^ : by tr in all terms in the series (B12). So 
tr Hl ln (y-j t ^ +iB^j = i £ dt G e (0) tr B (t) 

+1 f 1 dt f 1 dt' G e {t - t')G e (t' - t) trB(t)B{t') + ... (514) 
2 Jo Jo 

In the limit e — > +0 the first term in R.H.S. of (B14) is equal to 

- f 1 dttrB(t) 
2 jo 

(see (Bll)) whereas the sum of all others coincide with 

tr *' (-Jt + lB ) 

by virtue of (BIO'), (Bll'). This proves the validity of the equation (B2). 

To check the formula (B2), let us consider the simplest case of one 
dimensional harmonic oscillator. Partion function 

Z{(3) = trexp{-/3a t a} (515) 

is known exactly: 

OO 1 

Z{0) = £ e"^ = — — (515) 

n=0 1 e 

On the other hand, 

Z(P) = lim J pBc DzDz exp j jf * dt(-zz - Pze~ e *z J 

Putting in (B2) 5 = i(3, one gets: 
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(-1 (-I-*) (B17) 

The latter determinant we have already evaluated in the main text of 
the paper : 

det (--^ - pj = const sinh ^ (£18) 

(see eq. Q57|) for N = 1, r] = e = 0). 
Comparing (B16)-(B19), we get 

= const -^—s (£19) 
sinh | 



Evaluating the constant from the condition 



Z(oo) 



we reproduce the true answer (B15). We would like to stress that if we 
had put e = in (B16) before evaluating of the functional integral the 
result would have been wrong. 
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